CLNS 08/2019 
CLEO 08-03 



Two-Photon Widths of the Xcj States of Charmonium 

K. M. Ecklund, 1 W. Love, 2 V. Savinov, 2 A. Lopez, 3 H. Mendez, 3 J. Ramirez, 3 J. Y. Ge, 4 
D. H. Miller, 4 I. P. J. Shipsey, 4 B. Xin, 4 G. S. Adams, 5 M. Anderson, 5 J. P. Cummings, 5 
I. Danko, 5 D. Hu, 5 B. Moziak, 5 J. Napolitano, 5 Q. He, 6 J. Insler, 6 H. Muramatsu, 6 
C. S. Park, 6 E. H. Thorndike, 6 F. Yang, 6 M. Artuso, 7 S. Blusk, 7 S. Khalil, 7 J. Li, 7 
R. Mountain, 7 S. Nisar, 7 K. Randrianarivony, 7 N. Sultana, 7 T. Skwarnicki, 7 S. Stone, 7 

J. C. Wang, 7 L. M. Zhang, 7 G. Bonvicini, 8 D. Cinabro, 8 M. Dubrovin, 8 A. Lincoln, 8 
P. Naik, 9 J. Rademacker, 9 D. M. Asner, 10 K. W. Edwards, 10 J. Reed, 10 R. A. Briere, 11 
T. Ferguson, 11 G. Tatishvili, 11 H. Vogel, 11 M. E. Watkins, 11 J. L. Rosner, 12 
J. P. Alexander, 13 D. G. Cassel, 13 J. E. Duboscq, 13 R. Ehrlich, 13 L. Fields, 13 R. S. Galik, 13 

L. Gibbons, 13 R. Gray, 13 S. W. Gray, 13 D. L. Hartill, 13 D. Hertz, 13 J. M. Hunt, 13 
J. Kandaswamy, 13 D. L. Kreinick, 13 V. E. Kuznetsov, 13 J. Ledoux, 13 H. Mahlke-Kriiger, 13 
D. Mohapatra, 13 P. U. E. Onyisi, 13 J. R. Patterson, 13 D. Peterson, 13 D. Riley, 13 A. Ryd, 13 
A. J. Sadoff, 13 X. Shi, 13 S. Stroiney, 13 W. M. Sun, 13 T. Wilksen, 13 S. B. Athar, 14 
R. Patel, 14 J. Yelton, 14 P. Rubin, 15 B. I. Eisenstein, 16 I. Karliner, 16 S. Mehrabyan, 16 
N. Lowrey, 16 M. Selen, 16 E. J. White, 16 J. Wiss, 16 R. E. Mitchell, 17 M. R. Shepherd, 17 
D. Besson, 18 T. K. Pedlar, 19 D. Cronin-Hennessy, 20 K. Y. Gao, 20 J. Hietala, 20 Y. Kubota, 20 

T. Klein, 20 B. W. Lang, 20 R. Poling, 20 A. W. Scott, 20 P. Zweber, 20 S. Dobbs, 21 
Z. Metreveli, 21 K. K. Seth, 21 A. Tomaradze, 21 J. Libby, 22 A. Powell, 22 and G. Wilkinson 22 

(CLEO Collaboration) 

1 State University of New York at Buffalo, Buffalo, New York 14260, USA 
2 University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA 
3 University of Puerto Rico, Mayaguez, Puerto Rico 00681 
^ Purdue University, West Lafayette, Indiana 47907, USA 
5 Rensselaer Polytechnic Institute, Troy, New York 12180, USA 
6 University of Rochester, Rochester, New York 14627, USA 
7 Syracuse University, Syracuse, New York 13244; USA 
8 Wayne State University, Detroit, Michigan 48202, USA 

9 University of Bristol, Bristol BS8 1TL, UK 
10 Carleton University, Ottawa, Ontario, Canada K1S 5B6 
11 Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA 
12 Enrico Fermi Institute, University of Chicago, Chicago, Illinois 60637, USA 
ls Cornell University, Ithaca, New York 14853, USA 
1J > University of Florida, Gainesville, Florida 32611, USA 
15 George Mason University, Fairfax, Virginia 22030, USA 
16 University of Illinois, Urbana- Champaign, Illinois 61801, USA 
17 Indiana University, Bloomington, Indiana 47405, USA 
18 University of Kansas, Lawrence, Kansas 66045, USA 
19 Luther College, Decorah, Iowa 52101, USA 
20 University of Minnesota, Minneapolis, Minnesota 55455, USA 
21 Northwestern University, Evanston, Illinois 60208, USA 



1 



22 University of Oxford, Oxford 0X1 3RH, UK 
(Dated: March 19, 2008) 

Abstract 

Using a data sample of 24.5 million ip(2S) the reactions ip(2S) — ► 7Xcj> Xcj — * 77 have been 
studied for the first time to determine the two-photon widths of the Xcj states of charmonium 
in their decay into two photons. The measured quantities are B(ip(2S) — > 7Xco) x £>(XcO ~^ 
77 )=(2.22±0.32±0.10)xl0- 5 , and B(ip(2S) -► 7x c2 ) x B( Xc 2 77)=(2.70±0.28±0.15) x 10~ 5 . 
Using values for B(ip(2S) — > 7XcO,c2) and r(x c o,c2) from the literature the two-photon widths 
are derived to be r 77 (x c0 )=(2.53±0.37±0.26) keV, r 77 (x c2 )=(0.60±0.06±0.06) keV, and K = 
r 77 (xc2)/r 77 (xco)=0.237±0.043±0.034. The importance of the measurement of 1Z is emphasized. 
For the forbidden transition, Xci — ► 77, an upper limit of r 77 (x c i) < 0.03 keV is established. 
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Charmonium spectroscopy has provided some of the most detailed information about 
the quark-antiquark interaction in Quantum Chromodynamics (QCD). The most practical 
and convenient realization of QCD for onium spectroscopy is in terms of perturbative QCD 
(pQCD), modeled after Quantum Electrodynamics (QED). Two-photon decays of charmo- 
nium states Xcj^Pj) offer the closest parallel between QED and QCD, being completely 
analogous to the decays of the corresponding triplet states of positronium. Of course, the 
masses of the quarks and the wave functions of the Xc states differ from those of positronium, 
but even these cancel out in the ratio of the two-photon decays, so that for both positronium 
and charmonium TZ = L( 3 P 2 —>■ 77)/L( 3 Po ~~ > 77)=4/15~0.27 The departure from this 
simple lowest order prediction can arise due to strong radiative corrections and relativistic 
effects, and the measurement of TZ provides a unique insight into these effects. Two-photon 
decay of the spin one Xci state is forbidden by the Landau- Yang theorem 0]. There are 
numerous theoretical potential model predictions of L 77 (x c0iC2 ) available in the literature, 
with some employing relativistic and/or radiative corrections. As shown in Table [U the 
predictions vary over a wide range. This underscores the importance of measuring these 
quantities with precision. 

TABLE I: Potential model predictions for two-photon widths of Xc2 and XcO and the ratio TZ 
derived from them. 



Reference 


r 77 (x c2 ) (ev) r 


77(Xco) (eV) 


K 


Barbieri [3] 


930 


3500 


0.27 


Godfrey [4] 


459 


1290 


0.36 


Barnes [5] 


560 


1560 


0.36 


Bodwin [6] 


820±230 


6700±2800 


12+ - 15 

u - iz -0.06 


Gupta [7J 


570 


6380 


0.09 


Miinz [8] 


440±140 


1390±160 


S2+ - 16 


Huang [9] 


490±150 


3720±1100 


u - io -0.06 


Ebert [10] 


500 


2900 


0.17 


Schuler [11] 


280 


2500 


0.11 



Most of the existing measurements of L 77 (x c0 ) and L 77 (x c2 ) are based on formation of 
Xcj in two-photon fusion. The only existing measurements based on the decay of Xcj into 



two photons are from the Fermilab E760 /E835 experiments [12, ll3j, |14] with Xcj formation 
in pp annihilation. We report here results for L 77 (x c j) measured in the decay of Xcj into 
two photons. For these measurements we use the reactions 

ip(2S) -> 7!X c j, Xcj -> 727s, (1) 

which have not been studied before. Since L 77 (x c o) and L 77 (x c2 ) are obtained from the 
same measurement, we also obtain TZ with a good control of systematic errors. Few such 
simultaneous measurements have been reported in the literature. 

A data sample of 24.5 million ip(2S) obtained in 48pb~ 1 e + e~ annihilations at the CESR 
electron-positron collider was used. The reaction products were detected and identified using 
the CLEO-c detector. 



The CLEO-c detector [15j, which has a cylindrical geometry, consists of a Csl electromag- 



netic calorimeter, an inner vertex drift chamber, a central drift chamber, and a ring-imaging 
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Cherenkov (RICH) detector, inside a superconducting solenoid magnet providing a 1.0 T 
magnetic field. For the present measurements the most important component of the detector 
is the Csl calorimeter which has an acceptance of 93% of 47r and photon energy resolutions 
of 2.2% at £ 7 =1 GeV, and 5% at 100 MeV. 

The event selection for the final state required three photon showers, each with E 7 > 70 
MeV and angle 9 with respect e + beam direction with |cos#| < 0.75, and no charged 
particles. An energy-momentum conservation constrained 4C-fit was performed and events 
with x 2 /d.o.f. < 6 were accepted. To prevent overlap of the lowest energy photon 71 with 
the high energy photons 72,3, events were rejected if cos#' > 0.98, where 9' is the laboratory 
angle between 71 and either 72 or 73. 

Data were analyzed in two equivalent ways, by constructing the energy spectrum of £^(71) 
and the invariant mass spectrum of M(7 2 73). Consistent results were obtained. Fig. 1 shows 
the -E"(7i) spectrum. The enhancements due to the excitation of Xco and Xc2 over substantial 
backgrounds are clearly observed. 

In order to analyze these spectra we need to determine the shapes of the background 
and the resonance peaks. For determining peak shapes and efficiencies fifty thousand signal 
Monte Carlo (MC) events were generated for Xco and Xc2 each, with masses and widths as 
given by PDG 07 [la ]. The radiative transition if>(2S) — > 7iXco is, of course, pure El, and 
there is strong experimental evidence that the radiative transition ip(2S) — > 7iXc2 is also 
almost pure El 17, lj^. Further, 7973 in the decay Xc2 — ► 7273 are expected to be produced 



with pure helicity two amplitudes [5j. With these assumptions the angular distributions are 
predicted to be [19 



Xco ■ dN/ dcos&i = 1 + cos 2 Q u (2) 
Xc2 ■ d 3 N/(d cos ® t d cos 9 2 d0 2 ) = 9 sin 2 9i sin 2 29 2 
+ (1 + cos 2 00 [(3 cos 2 6 2 - l) 2 + 9 sin 4 Q 2 ] 
+3sin29i sin26 2 [3cos 2 9 2 — 1 — 3 sin 2 2 ] cos 02 
+6sin 2 eisin 2 e2(3cos 2 2 - l)cos20 2 , (3) 
Xc2 ■ dN/ dcosBi = 1 - (1/13) cos 2 0i. (4) 

Here 0i is the angle between 71 and the e + beam direction in the ip(2S) frame, and ©2 and 
<p2 are the polar and azimuthal angles of the 7273 axis in the rest frame of Xco,c2- These 
angular distributions were assumed in the MC simulations. 

The energy resolutions determined by the MC simulations were cr(-E 7l ) = (8.2±0.1) MeV 
for Xco and <r(i? 7l ) = (6.3±0.1) MeV for Xc2- The overall efficiencies determined from these 
MC samples were e(x c o) = (39.1±0.5)% and e(x c2 ) = (50.7±0.7)%. The difference between 
e (Xco) and e(x C 2) arises primarily from the cos 61 distributions (Eqs. (2) and (4)). 

Because the background in our spectrum is large, it is important to determine its shape 
accurately. For this purpose the distributions of £(71) were examined in the 21 pb _1 of 
o&-ip(2S) data taken at ^=3671 MeV, as well as the 280 pb -1 of large statistics -^(3770) 
data taken at ^/i=3772 MeV. As shown in Fig. 2, it is found that the o&-ip(2S) data are 
in excellent agreement with the high statistics ^(3770) data, in which transitions to either 
Xco or Xc2 resonances were expected to yield <2 events [20]. The £(71) distribution for the 
^(3770) data was fitted with a polynomial, and used as the shape of the background in the 
ip(2S) data shown in Fig. 1. 

The final fit to the E(ji) spectrum obtained with fixed Xco and Xc2 masses and intrinsic 
widths ] 1() is shown in Fig. 1. The x 2 /d-o.f. = 41/61 for the fit. It was found that 
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FIG. 1: Fitted spectrum for E(ji). The expected positions of E(ji) from XcO, Xci, Xc2 are marked 
with arrows. 
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FIG. 2: Background spectrum for E{^{). The points are from the off-ip(2S) data. The curve is 
from fit to the "0(3770) data. 

iV(x c o)=212±31 and iV(xc2)=335±35. The product branching ratios were determined as 
N{Xc.j)/Hxcj) x N(^(2S)} with the results 

Bty(2S) - 7Xco) x B( X co - 77) 

= (2.22 ± 0.32(stat)) x lfT 5 , . , 

BmS) -> 7Xc 2 ) x B( X c2 - 77) [ } 

= (2.70 ± 0.28(stat)) x 10" 5 . 

Various possible sources of systematic errors in our results were investigated. The number 
of ip{2S) produced is determined using the background-subtracted and efficiency-corrected 
yield of hadronic events following the procedure described in detail in [2l[. The background 
is estimated using the oS-i/j(2S) data. The efficiency is estimated by a MC simulation of 
generic ip{2S) decays. The systematic uncertainty is determined by varying the hadronic 
event selection and online trigger criteria by large amounts in both data and MC. It is 
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found that while the MC determined efficiency changes from 65% to 91% the efficiency- 
corrected yield changes by no more than 2%, which we include as a systematic error. The 
neutral trigger efficiency is uncertain by 0.2%. The uncertainty in our MC determination 
of absolute efficiency for three-photon detection was estimated as 3 x 0.4%=1.2% [221. The 
systematic error due to the simulation of the event selection criteria (x 2 /d.o.f. distribution 
for 4C fit, acceptance variation, and shower overlap cut) was determined by varying the 
cuts. Similarly, systematic uncertainties due to our choice of the background and signal 
shapes were estimated by using a free parameter polynomial background shape and a free 
parameter Crystal Ball line shape [2lJ convoluted with appropriate Breit-Wigner resonance 
shapes for the peaks. The extreme changes in the resonance yields obtained with these 
changes were taken as measures of systematic errors. We have assumed pure helicity two 
decay of Xc2- in a relativistic calculation Barnes 0] predicts the helicity zero component 
to be only 0.5%. To be very conservative, we have determined the change in our result for 
Xc2 by including a helicity zero component of 8%, which is the experimental upper limit 
established by CELLO [23] for the two photon decay of the 2 ++ light quark state 02(1320). 
All individual systematic errors are listed in Table [TT1 The sums of the systematic errors, 
added in quadrature are ±4.5% for Xco and ±5.7% for Xc2- 



TABLE II: Estimates of systematic uncertainties. Asterisks denote the systematic sources common 
to both Xco and Xc2- 



Source of Systematic Uncertainty Xco Xc2 



Number of ip(2S)* 


2.C 


% 


2.0 


% 


Neutral Trigger Efficiency* 


0.2 


% 


0.2 


% 


Photon Detection Efficiency * 


1.2 


% 


1.2 


% 


Event Selection Simulation 


2.C 


% 


2.0 


% 


Resonance Fitting 


3.2 


% 


4.6 


% 


Helicity 2 Angular Distribution 






1.3 


% 


Sum in quadrature 


A.l 


% 


5.7% 



Our final results for the measured quantities, B(ip(2S) — > 7Xco,c2) x B(xco,c2 — ► 77) are 
presented in Table IIHI We use the PDG 07 average results, 

(9.2 ±0.4) x 10~ 2 , 
(10.5 ±0.9) MeV, 
(8.8 ±0.5) x 10~ 2 , 

(1.95 ±0.13) MeV, (6) 

to derive B(xco,c2 —> 77), ^-yyiXco^), and 71. These are also listed in Table fTTTl 

By requiring an additional resonance in the spectrum of Fig. 1 corresponding to Xci{ 3 Pi), 
whose two-photon decay is forbidden by the Landau- Yang theorem [2|, we obtain the 90% 
confidence limit B(x c i —>■ 77) < 3.5 x 10 -5 , which is nearly two orders of magnitude lower 
than the present limit quoted in PDG 07 [16j. It corresponds to the 90% confidence limit 
r 77 (x c i) < 0.03 keV. 

Our final results are compared to those of previous measurements in Table IIVI As men- 
tioned earlier, most of the results for r 77 (x c j) in Table HVl are from measurements of the 



r(x c o) = 

B(^(2S) - 7X C2 ) = 

r(x c2 ) = 
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TABLE III: Results of the present measurements. First error is statistical, second is systematic, 
and third is due to the PDG parameters used. The common systematic errors have been removed 
in calculating TZ. Bi = B(ip(2S) -> 7XcO,c2), B 2 = B(xc0,c2 -> 77), T 77 = r 77 (x c o, C 2 -> 77)- 

Quantity XcO Xc2 

B 1 xB 2 x 10 5 2.22±0.32±0.10 2.70±0.28±0.15 
B 2 x 10 4 2.41±0.35±0.11±0.10 3.06±0.32±0.17±0.17 
T 77 (keV) 2.53±0.37±0.11±0.24 0.60±0.06±0.03±0.05 
K 0.237±0.043(stat)±0.015(syst)±0.031(PDG) 



formation of Xcj in two-photon fusion. The results listed in Table [IV] for r 77 (Xcj) have 
been updated by using the current PDG [l6[ values for the branching fractions and widths 
required for evaluating r 77 (x c j) from the directly measured quantities. 

From Table HV] we notice that although the individual results for r 77 (Xco) and r 77 (x c2 ) 
show large variations, their average agrees with our results. 

We note in Table [IV] that there are few simultaneous measurements of r 77 (Xco) and 
r~77 {Xcst)i an d even in those few cases TZ is seldom calculated. To put our result for TZ 
in perspective, we have calculated TZ for all simultaneous measurements in Table IIVI A 
weighted average of these is 7?.=0.20±0.03 which is in good agreement with our result, 
ft=0.24±0.06. 

The theoretical pQCD prediction for TZ is 32l 



TZ th = (4/15) [1 - 1.76a J , since 
r 77 ( Xc2 ) = 4(|*'(0)| 2 aL/m 4 ) x [1 - 1.70aJ , 
r 77 (Xco) = 15(|^'(0)| 2 aLMc) x [l + 0.06aj. 

where the quantities in the square brackets are the first order radiative correction factors. 



The radiative correction factor for r 77 (x c2 ) (for a s ~0.32 [16|) is nearly a factor of two, 
which strongly suggests possible problems with the radiative corrections. Unfortunately, a 
measurement of r 77 (x c2 ) alone cannot provide further insight into the problem because the 
charm quark mass m c and derivative of the wave function at origin \&'(0) are not known. 
However, since both unknowns cancel in the ratio TZ, a measurement of TZ can do so, as 



noted, for example, by Voloshin [33J. For a s =0.32, the predicted value, which only depends 
on radiative corrections, is TZ t h=0.l2. Our experimental result, 7£ exp =0.24±0.06, together 
with the other determinations of TZ in Table IIVI leads to the average (7£ exp )=0.20±0.02. 
This result provides experimental confirmation of the inadequacy of the present first-order 
radiative corrections, which have been often used to make theoretical predictions of r 77 (x c j) 
and experimental derivations of a s . 

The above experimental results for TZ emphasize the need for calculations of radiative 



corrections to higher orders. Alternatively, as noted by Buchmiiller |34j], a different choice 
of the renormalization scheme and renormalization scale should be considered in order to 
arrive at a more convergent way of specifying the radiative corrections. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the A. P. Sloan Foundation, 
the National Science Foundation, the U.S. Department of Energy, the Natural Sciences and 
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TABLE IV: Compilation of experimental results for two-photon partial widths of XcO an d Xc2- 



Experiment [Ref.] 


Measured 


r 77 ( Xc0 ) keV* 


r 77 ( Xc2 ) keV* 


K 


E760(1993) [12] 


£(pp -> Xc2) 


X 13-y-y 




0.47±0.12±0.07 




H/OOO^ZUUUJ 


£(pp -> Xcj) 


X 13 sysy 


o ni+1 n^-i-n 94 
z . u i in i . uoin u . Z4 


n QQ-i-n c\7-\-c\ 


n 90+0 1 1 -i-n 


TTSQ^/^nn/A ri/il 

-CjOO0^ZUU4J 1 141 


£>(pp — ► Xco) 


X 13-yy 

I I 


O.O±U.D±U.0 






OPAL (1998) [24] 


r 77 x fi(x C 2 - 


> 7 J/^) 


- 


1.19±0.32±0.26 


- 


TQ/1 QQQ\ rOC;l 


r 77 x fi(x C 2 - 


> 7-W 




u.oyztu.z i it u. ri 






r 77 x B(xc2 - 


> 7 J/V) 




U. 1 4t±U.Zl±U.lo 




UL_hj(J(2UUlJ 127 1 


T~~ X B (X f ] ~ 


* 7J/-0) 

/ It/ 


o.U9±U.54±U.44 


n £i _i_n 1 /)_i_n nn 

u.oi±u.i4±u.uy 


U.l /±U.Ud±U.U4 


CLEO(2006) [28] 


r 77 x B(xc2 - 


> 7 J/V) 




0.55±0.06±0.05 




Belle(2002) [29] 


r 77 x B(xc2 - 


> 7 J/V) 




0.56±0.05±0.05 




Belle(2007) [30] 


r 77 x B(xcj -» 


• *M 


2.53±0.23±0.40 


0.46±0.08±0.09 


0.18±0.03±0.04 


Belle(2007) [31]** 


r 77 x b(xc.j 


-> 4vr) 


1.84±0.15±0.27 


0.40±0.04±0.07 


0.22±0.03±0.05 




r 77 x B(xcj - 


► 2vr2if) 


2.07±0.20±0.40 


0.44±0.04±0.16 


0.21±0.03±0.09 




r 77 x e(x c j - 


-> 4if) 


2.88±0.47±0.53 


0.62±0.12±0.12 


0.21±0.05±0.06 


This measurement 


B^(2S) - 7X, 


; j) X £? 77 


2.53±0.37±0.26 


0.60±0.06±0.06 


0.24±0.04±0.03 


Averages (wei; 


2;hted by total errors) 


2.31±0.10±0.12 


0.51±0.02±0.02 


0.20±0.01±0.02 



* The first error is statistical. The second error is systematic error combined in quadrature with 
the error in the branching fractions and widths used. The results from the literature have been 
reevaluated by using the current PDG values for branching fractions and total widths. For these 
results the errors in 1Z have been evaluated without taking into account possible correlations in 
the systematic errors in r 77 (x c o) and r 77 (x C 2)- 

** The Belle publication gives only the product branching fractions T 77 x B(xc0,c2 — ► hadrons). 
We have calculated T 77 and 1Z by using the PDG 07 [16| values of branching fractions for the 
individual decays. 

Engineering Research Council of Canada, and the U.K. Science and Technology Facilities 
Council. 



[1] See, for example, V.A. Novikov et al, Phys. Rep. 41, 1 (1978). 

[2] L. Landau, Phys. Abstracts A 52, 125 (1949); C.N. Yang, Phys. Rev. 77, 242 (1950). 

[3] R. Barbieri, R. Gatto and R. Kogerler, Phys. Lett. B 60, 183 (1976). 

[4] S. Godfrey and N. Isgur, Phys. Rev. D 32, 189 (1985), as quoted in Ref. [7]. 

[5] T. Barnes, Proceedings of the IX International Workshop on Photon-Photon Collisions, edited 

by D.O. Caldwell and H.P. Paar (World Scientific, Singapore, 1992), p. 263. 
[6] G. Bodwin, E. Braaten and G.P. Lepage, Phys. Rev. D 46, R1914 (1992). 



8 



[7] S.N. Gupta, J.M. Johnson and W.W. Repko, Phys. Rev. D 54, 2075 (1996). 
[8] C.R. Miinz, Nucl. Phys. A 609, 364 (1996). 

[9] H.-W. Huang and K.-Ta Chao, Phys. Rev. D 54, 6850 (1996); errata Phys. Rev. D 56, 1821 
(1997). 

[10] D. Ebert, R.N. Faustov and V.O. Galkin, Mod. Phys. Lett. A 18, 601 (2003). 
[11] G.A. Schuler, F.A. Berends and R. van Gulik, Nucl. Phys. B 523, 423 (1998). 
[12] T.A. Armstrong et al. (E760 Collaboration), Phys. Rev. Lett. 70, 2988 (1993). 
[13] M. Ambrogiani et al. (E835 Collaboration), Phys. Rev. D 62, 052002 (2000). 
[14] M. Andreotti et al. (E835 Collaboration), Phys. Lett. B 584, 16 (2004). 

[15] G. Viehhauser, Nucl. Instrum. Methods Phys. Res., Sect. A 462, 146 (2001); D. Peterson et 

al, Nucl. Instrum. Methods Phys. Res., Sect A 478, 142 (2002). 
[16] W.-M. Yao et al, Journal of Physics G 33, 1 (2006) and 2007 partial update for the 2008 

edition available on the PDG WWW (URL: http://pdg.lbl.govD . 
[17] M. Oreglia et al. (Crystal Ball Collaboration), Phys. Rev. D 25, 2259 (1982). 
[18] M. Ambrogiani et al. (E835 Collaboration), Phys. Rev. D 65, 052002 (2002). 
[19] PK. Kabir and A.J.G. Hey, Phys. Rev. D 13, 3161 (1976). 

[20] T.E. Coan et al. (CLEO Collaboration), Phys. Rev. Lett. 96, 182002 (2006); R.A. Briere et 

al. (CLEO Collaboration), Phys. Rev. D 74, 031106 (2006). 
[21] S.B. Athar et al. (CLEO Collaboration), Phys. Rev. D 70, 112002 (2004). 
[22] N.E. Adam et al. (CLEO Collaboration), Phys. Rev. Lett. 94, 232002 (2005). 
[23] H.J. Behrend et al. (CELLO Collaboration), Z. Phys. C 46, 583 (1990). 
[24] K. Ackerstolt et al. (OPAL Collaboration), Phys. Lett. B 439, 197 (1998). 
[25] M. Acciarri et al. (L3 Collaboration), Phys. Lett. B 453, 73 (1999). 
[26] J. Dominick et al. (CLEO Collaboration), Phys. Rev. D 50, 4265 (1994). 
[27] B. Eisenstein et al. (CLEO Collaboration), Phys. Rev. Lett. 87, 061801 (2001). 
[28] S. Dobbs et al. (CLEO Collaboration), Phys. Rev. D 73, 071101 (2006). 
[29] K. Abe et al. (Belle Collaboration), Phys. Lett. B 540, 33 (2002). 
[30] W.T. Chen et al. (Belle Collaboration), Phys. Lett. B 651, 15 (2007). 
[31] S. Uehara et al. (Belle Collaboration), Eur. Phys. J. C 53, 1 (2008). 

[32] R. Barbieri et al, Phys. Lett. B 95, 93 (1980); Nucl. Phys. B 192, 61 (1981). For a review 

see W. Kwong et al, Phys. Rev. D 37, 3210 (1988). 
[33] M.B. Voloshin, arXiv:hep-ph/0711.4556. 

[34] W. Buchmuller, in "Quarkonia", in Current Physics Sources and Comments", vol. 9, North 
Holland Publishers, 1992. 



9 



